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Abstract 



We compute quark structure functions and the intrinsic Fock space distribution 
of sea quarks in a hadron wavefunction at small x. The computation is performed 
in an effective theory at small x where the gluon field is treated classically. At 
Q 2 large compared to an intrinsic scale associated with the density of gluons /x 2 , 
large compared to the QCD scale Aq CD , and large compared to the quark mass 
squared M 2 , the Fock space distribution of quarks is identical to the distribution 
function measured in deep inelastic scattering. For Q 2 < M 2 but Q 2 » fi 2 , 
the quark distribution is computed in terms of the gluon distribution function and 
explicit expressions are obtained. For Q 2 < [i 2 but Q 2 » A.q CD we obtain formal 
expressions for the quark distribution functions in terms of the glue. An evaluation 
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of these requires a renormalization group analysis of the gluon distribution function 
in the regime of high parton density. For light quarks at high Q 2 , the DGLAP 
flavor singlet evolution equations for the parton distributions are recovered. Explicit 
expressions are given for heavy quark structure functions at small x. 
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1 Introduction 



One of more interesting problems in perturbative QCD is the behaviour of structure 
functions at small values of Bjorken x. In deep inelastic scattering (DIS) for instance, 
for a fixed Q 2 >> Aq CD , the operator product expansion (OPE) eventually breaks 
down at sufficiently small x JjJ]. Therefore at asymptotic energies, the conventional 
approaches towards computing observables based on the linear DGLAP Q equations 
are no longer applicable and novel techniques are required. Even at current colider 
energies such as those of HERA where the conventional wisdom is that the DGLAP 
equations successfully describe the data, there is reason to believe that effects due to 
large logarithms in as log(l/x) (or large parton densities) are not small and we may 
be at the threshold of a novel region where non-linear corrections to the evolution 
equations are large ||, |j. One reason violations of DGLAP evolution have not been 
seen clearly thus far at HERA is the small phase volume for gluon emission ||. A 
straightforward way to further probe this region at current colider energies would 
be by using nuclear beams at HERA Q or in electron-proton collisions at the LHC 
collider where the phase volume for gluon emission is signficantly larger J?]]. 

In recent years, a non-OPE based effective field theory approach to small x 
physics has been developed by Lipatov and collaborators ||. Their initial efforts 
resulted in an equation known popularly as the BFKL equation ||, which sums 
the leading logarithms of aslog(l/x) in QCD. In marked contrast to the leading 
twist Altarelli-Parisi equations for instance, it sums all twist operators that contain 
the leading logarithms in x. The solutions to the BFKL equation predict a rapidly 
rising gluon density and there was much initial euphoria when the HI and ZEUS 
data at HERA showed rapidly rising parton densities However, it was shown 
since: 

• The rapid rise of the structure functions can arguably be accounted for by 
the next to leading order (NLO) DGLAP equations by appropriate choices of 
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initial parton densities [11, 12 1. 



• The next to leading logarithmic corrections to the BFKL equation computed in 



the above mentioned effective field theory approach are very large |13L 14, 



As a consequence, the theoretical situation is wide open and novel approaches need 
to be explored. 

An alternative effective field theory approach to QCD at small x was put forward 
in a series of papers p6[-p0|]. In the approach of Lipatov and collaborators, the 
fields of the effective theory are composite reggeons and pomerons. This approach 
is motivated by the reggeization of the gluon that occurs in the leading log result. 



Our approach based on Refs.[16]-[20| is instead a Wilson renormalization group 
approach where the fields are those of the fundamental theory but the form of 
the action at small x is obtained by integrating out modes at higher values of x. 
Integrating out the higher x modes results in a set of non-linear renormalization 



group equations [20]. If the parton densities are not too high, the renormalization 



group equations can be linearized and have been shown to agree identically with 



the leading log BFKL and small x DGLAP equations [19]. There is much effort 
underway to explore and make quantitative predictions for the non-linear regime 
beyond || H- 

In this paper, we apply the above effective action approach to study the fermionic 
degrees of freedom at small x. At small values of x, gluon degrees of freedom domi- 
nate and the fermionic degrees of freedom present are essentially the sea quarks that 
are radiatively generated from the glue |?4|] (and are therefore O(as) corrections). 
Nevertheless, the sea quark distributions are extremely important since they are 
directly measured in deep inelastic scattering experiments. In this paper, we will 
develop a formalism, in the context of the renormalization group approach, which 
relates structure functions at small x to the sea quark distributions, and therefore to 
the gluon distribution. We derive analytic expressions summing a particular class 
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of all twist operators which we argue give the dominant contribution at small x. 
At leading twist, in light cone gauge, these reduce to the well known simple rela- 
tion between the structure function F<i and the sea quark Fock space distribution 
function 25] 

i dN SKn 



F 2 {x,Q 2 ) = I dk 
Jo 



2'_ 

* dk 2 dx 



Above, x and Q 2 are the usual invariants in deep inelastic scattering. We show 
explicitly that for light quarks and at high Q 2 , we reproduce the Altarelli-Parisi 
evolution equations for the quark distributions at small x. 

There has been much interest in heavy quark distributions motivated partly by 
the significant contribution of heavy quarks to the structure functions at HERA |30| 1 . 
Until recently, heavy quarks were treated as infinitely massive for Q 2 equal to or 
less then the quark mass squared and massless below. There now exist approaches 
which study quark distributions in a unified manner for a range of Q 2 and quark 



masses (for a discussion and further references see Refs. [31, |3^]). A nice feature 



of our formalism is that heavy quark evolution is treated on the same footing as 
light quarks and specific predictions can be made for p^ iarm / bottom j p 2 within our 
formalism. These can also be related to the diffractive cross section at small x but 
that issue will not be addressed in this paper. This issue and the relation of our 
work to the above cited work and other recent works on heavy quark production at 



small x |34|] will be addressed at a later date. 

The results of our analysis are the following. In our theory of the gluon distri- 
bution functions, a dimensionful scale appears which measures the density of gluons 
per unit area, 

a dy 

where a is the hadronic cross section of interest. Here y = yo — ln(l/x), yo is 
an arbitrarily chosen constant and x is Bjorken x. When this parameter satisfies 



fi » Aqcd the gluon dynamics, while nonperturbative, is both weakly coupled and 
semiclassical. We shall always assume that at sufficiently small x this is satisfied. 

At small x, the gluon field, being bosonic, has to be treated non-perturbatively. 
This is analogous to the strong field limit used in Coulomb problems. Fermions, 
on the other hand, do not develop a large expectation value and may be treated 
perturbatively. To lowest order in uqcd, the gluon distribution function is deter- 
mined by knowing the fermionic propagator in the classical gluon background field. 
In general, this propagator must be determined to all orders in the classical gluon 
field as the field is strong. This can be done due to the simple structure of the 
background field. 

At this point it is useful to distinguish between two different quantities which 
are often used interchangeably. One is the quark structure functions as measured 
in deep inelastic scattering and the second is the Fock space distribution of quarks. 
At high Q 2 which is usually the case considered in perturbative QCD, these two 
quantities are essentially identical. However, for massive quarks when Q 2 < M 2 
even if Q 2 » Aq CD , the two quantities differ. For Q 2 » fi 2 , regardless of mass, 
the quark Fock space distribution function and the quark structure functions may 
both be simply expressed as linear functions of the gluon distribution functions. 
Therefore, with knowledge of the gluon distribution function, one can compute the 
quark distribution function. 

For the case of massless quarks, when y? > Q 2 » Aq CD , we are still in the 
weak coupling limit. However, we must keep all orders in the gluon field. In this 
region, the integration over the gluon fields in our effective field theory cannot be 
directly performed as yet since it requires a full renormalization group analysis of 
the theory. In other words, the measure of integration for the high density regime in 
the effective theory has not yet been computed. Nevertheless, we obtain an explicit 
functional dependence on the gluon fields which must be integrated over with the 
right measure. 
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The power of the technique which we use to analyze this problem is that it 
does not rely on a high twist expansion. It uses only the weak coupling nature of 
the theory which must be true at small x if the gluon density is very high. We are 
therefore in a position to find non-trivial relations between these various parameters 
in a region where the weak coupling analysis is valid but where perturbation theory 
and leading order operator product expansion methods are not valid. 

This paper is organized as follows. In section 2 we write down and review an 
effective action for the small x modes in QCD. The action is imaginary and the 
modes are averaged over with a statistical weight 

exp [-F[p}\ , 

where F[p] is a functional over the color charge density p of the higher x modes. The 
functional F[p] obeys a non-linear renormalization group equation. Of particular 
interest in this paper is the saddle point solution of the effective action since the sea 
quark distributions are computed in the classical background field of this action. 
This section is a quick review of known results which are necessary to understand 
the remainder of the paper. 

An expression relating the electromagnetic current-current correlator to the 
fermion propagator in the classical background field is derived in section 3. We also 
discuss the light cone Fock space distributions and their relation to the structure 
functions in this section. 

In section 4, we solve the Dirac equation in the classical background field and 
obtain an explicit expression for the fermion Green's function in the classical back- 
ground field. 

The Green's function is used in section 5 to compute the sea quark distribution 
function and the leading twist contribution to the structure function F2(x,Q 2 ) at 
small x. The color averaging over the functional F[p] in the distribution function 
is compactly represented by a function ^(pt), where pt can be interpreted as the 
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intrinsic transverse momentum of the glue at high parton densities. It is shown 
explicitly that for light quark masses the flavor singlet Altarelli-Parisi equations at 
small x are recovered. 

The current -current correlator and the structure functions are computed explic- 
ity in section 6. As a check of our computation, it is shown that the leading twist 
results are recovered in the appropriate limit. The heavy quark structure functions 
are computed explicitly. The phenomenological implications and the connections to 
the recent literature on heavy quark production will not be addressed in this work 
but will be considered at a later date. 

Section 7 contains a summary of our results and a discussion of future work. 

The first of two appendices contains a discussion of our notation and conven- 
tions. In the second appendix we present an explicit form for y(pt) for the particular 
case of Gaussian color fluctuations. 

2 Effective Field Theory for Small x Partons in QCD: 
Review of Results 

We will discuss below an effective action for the wee parton modes in QCD. The 
action contains an imaginary piece which involves functional F[p] which satisfies a 
non-linear Wilson renormalization group equation. In the weak field limit of this 
renormalization group equation, the BFKL equation is recovered. In the double 
logarithmic region, the evolution equation is also equivalent to DGLAP Q. We next 
discuss the classical background field which is the saddle point solution of this 
action. It is this background field that the sea quarks couple to at small x and 
the properties of the background field will be relevant for the discussion in later 
sections. 

1 This was first noticed by Yuri Dokshitzer in his paper in Ref. M. 
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2.1 The Effective Action and the Wilson Renormalization Group 
at Small x 



In the infinite momentum frame P + — > oo, the effective action for the soft modes 
of the gluon field with longitudinal momenta k + « P + (or equivalently x = 
k + / P + « 1) can be written in light cone gauge A + = as 

Seff = - J ' d+x-G^G^ + ^- j ' d 2 x t dx- p a {x u x-)Tr{r a W- 00 , 00 [A-}{x--,x t )) 
+ i J d 2 x t dx-F[p a (x t ,x-)}. (2) 

Above, Gp V is the gluon field strength tensor, r a are the SU(N C ) matrices in the 
adjoint representation and W is the path ordered exponential in the x + direction 
in the adjoint representation of SU(N C ), 



W-oo )00 [A ]{x ,x t ) = Pexp 



ig / dx + A a (x , xt)r a 



(3) 



The above is the most general gauge invariant form [19] of the action that was 



proposed in Ref. [16| 



This is an effective action valid in a limited range of P + « A + where A + is an 
ultraviolet cutoff in the plus component of the momentum. The degrees of freedom 
at higher values of P + have been integrated out and their effect is to generate the 
second and third terms in the action. 

The first term in the above is the usual field strength piece of the QCD action 
and describes the dynamics of the wee partons at the small x values of interest. 
The second term in the above is the coupling of the wee partons to the hard color 
charges at higher rapidities, with x values corresponding to values of P + > A + . 
When expanded to first order in A~ this term gives the ordinary J ■ A coupling 
for classical fields. The higher order terms are needed to ensure a gauge invariant 
coupling of the fields to current. 
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In the infinite momentum frame, only the J + component of the current is large 
(the other components being suppressed by 1/P + ). The longer wavelength wee 
partons do not resolve the higher rapidity parton sources to within 1/P + and for 
all practical purposes, one may write 

p a (x t ,x-) ^ p a {x t )5{x-). (4) 

The last term in the effective action is imaginary. It can be thought of as a statistical 
weight resulting from integrating out the higher rapidity modes in the original QCD 
action. Expectation values of gluonic operators 0{A) are then defined as 

t r>( , 4 w _ J[dp}exp(-F[p}) f[dA]0(A)exp(iS[p,A]) 

{ ] ~ J[dp]ex V (-F[p])I[dA]ex V (iS[p,A]) ' 1 ' 

where S[p, A] corresponds to the first two terms in Eq. ^. The color averaging 
procedure for fermionic observables is discussed further in section 5.1. 
In Ref. [16] a Gaussian form for the action 

\ Sxt ^pl papa ' (6) 
was proposed, where p? was the average color charge squared per unit area of the 
sources at higher rapidities than is appropriate for our effective action, that is. For 
large nuclei A >> 1 it was shown that 

1 N n 



,, 2 1 
P = =- ^ 



AW/Gfmr 2 . (7) 



irR 2 2N c 

This result was independently confirmed in a model constructed in the nuclear rest 



frame [22]. If we include the contribution of gluons which have been integrated out 



by the renormalization group technique, one finds that [44| 



M^f^ + #^) (8) 



ttR 2 \2N c iV c 2 - 1 
Here N„ is the total number of quarks with x above the cutoff 



E f dx' qi (x') (9) 
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where the sum is over different flavors, spins, quarks and antiquarks. For gluons, 
we also have 



The value of irR 2 is well defined for a large nucleus. For a smaller hadron, we 
must take it to be a, the total cross section for hadronic interactions at an energy 
corresponding to the cutoff. This quantity will become better defined for a hadron 
in the renormalization group analysis. 

The above equation for /i 2 is subtle because, implicitly, on the right hand side, 
there is a dependence on fx through the structure functions themselves. This is the 
scale at which they must be evaluated. Calculating \x therefore involves solving an 
implicit equation. Note that because the gluon distribution function rises rapidly 
at small x, the value of fi grows as x decreases. 

The Gaussian form of the functional F[p] is reasonable when the color charges 
at higher rapidity are uncorrelated and are random sources of color charge. This 
is true for instance in a very large nucleus. It is also true if we study the Fock 
space distribution functions or deep inelastic structure functions at a transverse 
momentum scale which is larger than an intrinsic scale set by asH- In this equation 
as is evaluated at the scale fi. At smaller transverse momenta scales, one must do a 
complete renormalization group analysis to determine F[p\. This analysis is not yet 
complete, but should be feasible in the context of the weakly coupled field theory as 
long as the transverse momentum scale remains larger than that Aqcd- The color 
averaging procedure for the case of Gaussian fluctuations is discussed in appendix 



A final comment about the limit of applicability of the classical action above 
concerns limitations in the transverse momentum range. The action above is valid 
only when probing transverse momenta scales px < H- This includes the Gaussian 
region since as « 1. At higher transverse momenta, one must use DGLAP evolu- 




(10) 



B. 



11 



tion, with the structure functions as determined at lower values of Q 2 as boundary 



conditions [45]. This will be important for the case of heavy quarks as the transverse 
momentum scale there is very large. 

The above comments on the renormalization group analysis show the limitations 
of our analysis with respect to quarks. For transverse momentum scales pt >> a s ft, 
one can use a Gaussian source and all relevant quantities can be computed explicitly. 
At smaller scales, one can derive a formal expression, which, hopefully, will be 
directly computed in the near future. For heavy quarks, the Gaussian analysis 
should be adequate. 

To complete a review of the renormalization group, we briefly review the pro- 
cedure used to determine F[p]. It was shown that a Wilson renormalization group 
procedure [^] could be applied to derive a non-linear renormalization group equa- 
tion for F[p\. The procedure, briefly, is as follows. The gauge field is split as 

A»(x) = b»(x) + 6Afa) + a°(x), (11) 

where is the saddle point solution of Eq. [2] and corresponds to the hard 

modes above the longitudinal momentum scale A. The fluctuation field 5A^(x) 
contains the soft modes A /+ < k + < A + and a^{x) are soft fields with longitudinal 
momenta k + < A /+ . The cutoffs are chosen such that as log(A' + /A + ) << 1. Small 
fluctuations are performed about the saddle point solution b^(x) to the effective 
action at the scale A + , to obtain the effective action for the fields a^x) at the scale 
A /+ . The new charge density at this scale p' is given by p a ' = p a + 5p a , where 5p a 
can be expressed as the sum of linear and bi-linear terms in the fluctuation field 
6A«(x). 

To leading order in as, the effective action at the scale A /+ can be expressed in 
the same form as Eq. |2|, with the functional F[p'] satisfying a non-linear renormal- 
ization group equation [pll]-[2C]. In terms of the statistical weight Z = exp(— F[p\), 
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it can be expressed as 
dZ 



\ tTx ( z X^) ~ y~ (Zap) 
2 bp^dpv r bp^ 



(12) 



d\og{l/x) 

where a[p] and x[p] are respectively one and two point functions obtained by inte- 
grating over 5 A for fixed p. The one point function a includes the virtual corrections 
to F[p] while the two point function x includes the real contributions to F[p]. Both 
of these can be computed explicitly from the small fluctuations propagator in the 
classical background field. The propagator was first computed by fixing the residual 
gauge freedom to be 5A~ (x~ = 0) |l7j but a less restrictive gauge choice was later 



found which may be useful for computing a and X [20|. 

For weak fields, the free gluon propagator can be used to obtain the well known 
BFKL equation for the unintegrated gluon density, which is defined as 



dN 
dtf 



d xte 



-ik t x t 



(13) 



Performing the renormalization group procedure defined above to obtain the charge 
density p' = p + bp, one obtains 



< PP' >p~ < PP >p= "5 log(l/x) [2 < pa > p + < x > jC 



(14) 



The one and two point functions a and x respectively can be computed to linear 
order in the classical background field and the results are Jn 

9 2 N C 



<J a {h 



2(2vr) : 



■ P a (h 



d 2 pt 



PtiPt - h) 2 ' 



and 



X 



2g 2 N c 



d 2 PtP a (pt) P a (-Pt) 



kl 



(15) 



(16) 



(2vr)3 J — — ^ ^pKh-pt?' 
The above are respectively the virtual and real contributions to the change in the 
color change density after integrating out the modes A /+ < k + < A + . Substituting 
these into Eq. O, one obtains the well known BFKL equation 



dN 
dkfdx 



asNc 
2vr 2 



d p t7 2 



kf 



PtiPt 



h) 2 



dN 



dN 



dk\ dp 2 



(17) 
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Strenuous efforts are currently underway to compute a and x to all orders in 
the background field and thereby solve the full non-linear Wilson renormalization 
group equation for F[p] ^l], |23|| . 



2.2 The Classical Background Field at Small x 



The effective action in Eq. |2| has a remarkable saddle point solution [16, 18, |22j| . It 
is equivalent to solving the Yang-Mills equations 

= J U 5 U+ , (18) 

in the presence of the source J +,a = p a (xt, x~). Here we will allow the source to 
be smeared out in x~ as this is useful in the renormalization group analysis. It is 
also useful for intuitively understanding the nature of the field. One finds a solution 
where A = and 

A 1 = — Vd^ , (19) 

for % = 1, 2 is a pure gauge field which satisfies the equation 

dA i 

Di—=gp(y,x ± ). (20) 

Here Di is the covariant derivative di + VdiV* and y = yo + log (a; ~ /x^ ) is the 
space-time rapidity and yo is the space-time rapidity of the hard partons in the 
fragmentation region. At small x we will use the space-time and momentum space 
notions of rapidity interchangeably [45|. The momentum space rapidity is defined 



to be y = yo — ln(l/x) where x is Bjorken x. The solution of the above equation is 
A^ Xt ) = l( Pe 9J y V ^ y ' M v « ( Pe 9J y y ^ y > M (21) 
To compute the classical nuclear gluon distribution function, for instance, 

= (27j3 2|fc+l / d3xd3x ' elHx ~ X,) < A ?(x-^t)A«(x>-,x' t ) > p , (22) 
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one needs in general to average over the product of the classical fields at two space- 
time points with the weight F[p] as shown in Eq. [| or for the Gaussian measure 
with the weight in Eq. ||. In the latter case, exact analytical solutions are available 
for correlators in the classical background field. For the case of interest here, the 
fermion Green's function in the classical background field will depend on correlators 
of the form < V(xt)V'(yt) > p where the V's are SU(N C ) gauge transformation 
matrices defined above. In appendix B, we discuss in detail the computation of this 
correlator for the case of Gaussian fluctuations. 



3 The Current— Current Correlator at Small x 



In this section, we will derive a formal expression for the hadron tensor W pu at 
small x relating it to the fermion Green's function in the classical background field. 
We also derive a relation between the light cone quark distribution function and 
the fermion Green's function. To leading twist, the structure functions are simply 
related to the light cone quark distribution function. In general (for example, for 
heavy quark distributions) this is not true. Nevertheless, a simple relationship may 
be found between the gluon distribution functions and that of the quarks. This is 
because the quarks distribution functions is given by an integral over a propagator 
in the classical gluon background field described above. 

3.1 Derivation 



In deep inelastic electroproduction, the hadron tensor can be expressed in terms of 
the forward Compton scattering amplitude T^ u by the relation pq ] 

W^(q 2 ,P-q) = 2DiscT^(q 2 ,P-q) = — Im f d 4 x exp{iq ■ x) 

2ir J 

x < P\T(JV( x )J v (0))\P >, (23) 
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where "T" denotes time ordered product, = ip^ip is the hadron electromagnetic 
current and "Disc" denotes the discontinuity of T^ v along its branch cuts in the 
variable P ■ q. Also, q 2 — ► oo is the momentum transfer squared of the virtual 
photon ^] and P is the momentum of the target. In the infinite momentum frame, 
P + — > oo is the only large component of the momentum. The fermion state above 
is used in the expectation value for the current operators is normalized as < P \ 
P' >= (27r) 3 £^/m 5^ (P — P') where m is the mass of the target hadron. This 
definition of W^ v and normalization of the state is traditional, and we will abide by 
these conventions in spite of the awkward factors of m. We will see in the end that 
all factors of m cancel from the definition of quantities of physical interest. (The 
normalization we will use in this paper for quark and lepton states will have E/m 
replaced by 2P + .) 

Let us first describe the computation of < P | T(J^{x)J v (G)) \ P >. In our 
computation, we have an external source corresponding to the particle whose state 
vector is denoted by | P >. Our source is located at some fixed position. We must 
therefore consider the generalization of W^ v for such a source which has a position 
dependence. Note that for a given source, we also have a lack of translational 
invariance in the transverse direction. Transverse translational invariance is restored 
after integration over the source. There is no dependence of our source on x + . 
Therefore the relevant variable is x~ . 

We now argue that the relevant definition of W^ u is 

W^iq^P-q) = —a—Im 
Z7T m 

To see this let us first verify that we can write in this form for the conventional 
definition valid for plane wave states \P >. Notice that we can define < O >=< 
P | O | P > / < P \ P > where O is any operator. As mentioned above, the 
2 Note that in our metric convention, a space-like photon has q 2 = Q 2 > 0. 



J d 4 xdX~e iq - x < T (J»{X- + x/2)J u (X~ - x/2)) > . 
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expectation value < P \ P >= (2n) 3 E/m <5 (3) (0) = {2-nfEjm V. Here we shall 
take the spatial volume V to be a times an integral over the longitudinal extent of 
the state. Using these conventions, we see that we reproduce the above definition 
of W^. 

This definition corresponds to treating the variable X~ as a center of mass 
coordinate and X clS &■ relative longitudinal position. For a translationally invariant 
state, this would give the longitudinal dimension of the system. The definition is 
Lorentz covariant as will be shown explicitly in section 6. The integration over X~ 
is required since we must include all of the contributions from quarks at all X~ to 
the distribution function. In our external source language, the variable P + can be 
taken to be the longitudinal momentum corresponding to the fragmentation region. 
In the end all of the P + (and m) dependence will disappear upon taking the infinite 
momentum limit. To check this definition later, we shall show that it reproduces 
the conventional results in the high q 2 region. 

The expectation value is straightforward to compute in the limit where the 
gluon field is treated as a classical background field. If we write 

< T(J"(x)r(y)) >=< T (?(s) 7 ^0«#(l/)7 , V(y)) > , (25) 

then when the background field is classical, and one ignores quantum corrections 
arising from either loops of fermions or loops of gluons, (a good approximation in 
the weak coupling limit of high parton densities) , we obtain 

< T(J»(x)r(y)) >= Tr(^S A (x))TT(YS A (y)) + Tr(-fS A (x, y)-fS A (y, x)) . (26) 

In this expression, S A (x,y) is the Green's function for the fermion field in the 
external field A 

S A (x, y) = -i< ^{x)^(y) > A (27) 
for fixed A (before averaging over A). 
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The first term on the right hand side of Eqn. ^ is a tadpole contribution which 
does not involve a non-zero imaginary part. It therefore does not contribute to 
W^. We find then that 

W^(q 2 ,p-q) = —a—Im[dX-d A xe iqx (Tr(j fl S A (X~+x/2,X--x/2) 

Z7T 771 J V 

x -fS A (X--x/2,X-+x/2))). (28) 

3.2 The Light Cone Fock Distribution Function and Structure Func- 
tions 



The expression we derived above for is entirely general and makes no reference 
to the operator product expansion. In particular, it is relevant at the small x 
values and moderate q 2 where the operator product expansion is not reliable Q . At 
sufficiently high q 2 though (and for massless quarks) it should agree with the usual 
leading twist computation of the structure functions. The fact that we do not have 
a valid operator product expansion forces us to distinguish between two quantities 
which are identical in leading twist. The first is the Fock space distribution of 
partons within a hadron. The second are the parton structure functions which 
are measured in deep inelastic scattering. In our analysis, at high values of q 2 , 
these expressions are identical. At smaller values, say those values typical of the 
intrinsic transverse momentum scale cts/i, they are no longer the same and must be 
differentiated between. 

We will derive below an expression for the light cone quark Fock distribution in 



terms of the propagator in light cone quantization [4C, 41]. The quark Fock space 
distribution is then simply related to the structure function Fi for q 2 » a^fi 2 . In 
light cone quantization, only the two component spinor projection tp + is dynamical. 
(Note: notation and conventions are discussed in appendix A.) The other two spinor 
components ip- (recall that ip = tp_ + are defined via the light cone constraint 



relation defined below in Eq. 41. The dynamical fermions can then be written in 
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terms of creation and annihilation operators as 



d 3 k 



k+>0 



2l/4(27r)3 E [e tk - x w(s)b s (k)+e-^w(s)di(k) 



(29) 



Above b s (k) is a quark destruction operator and destroys a quark with momentum 
k while d\(k) is an anti-quark creation operator and creates an anti-quark with 
momentum k. Also above the unit spinors w(s) are defined as 



«4> 



1 







1 



(30) 



Note that since 



S(k t ,k + ;x + ;+-) 
S(k u k + ;x + ;- 1 -) 



b(-k t ,-k + ;x + ;--) 
b(-k t ,-k + ;x + ;+-) 



one can show that 



w(s)b s (k) = 2 1 ' i J d 3 xe- ik - x ip+, s (x). 



(31) 



The light cone Fock distribution function is defined in terms of the creation and 
annihilation operators as 



btb 



s u s 



(32) 



We have assumed above that the sea is symmetric between quarks and anti-quarks. 
Combining the two equations above, we get 
dN 2^2 



d 3 k (2vr) 3 
Using the light cone identity 



d 3 xd 3 ye ik < x - y ^l(x)ip + (y) 



Tr [7+^)^(2/)] = >/2 V+(aOVi(y) , 



(33) 



(34) 
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we obtain the following expression for the sea quark distribution function 
dN 2i 



d 3 x d 3 y e -*M Tr [ 7 +S A (x, y)] (35) 



d?k (2tt) 3 

where the fermion propagator S A {x,y) is the light cone time ordered product 
Sa(x,u) = —i < T(ip(x)ifj(y)) > in the background field A^. In our effective action 
approach, as discussed in sections 2 and 3.1, we can replace Sa(x,u) — ► Sa c1 (x,u) 
to obtain the sea quark distribution in the classical background field. 



In a nice pedagogical paper, (see Ref. [25] and references within), Jaffe has 
shown that the Fock space distribution function can be simply related to the leading 
twist structure function F% by the relation 

F 2 (x,Q 2 ) = j o dk 2 t x^-. (36) 

Actually, Jaffe's expression is defined as the sum of the quark and anti-quark dis- 
tributions. At small x, these are identical and the resulting factor of 2 is already 
included in our definition of the light cone quark distribution function. In the follow- 
ing section, we will compute the fermion Green's function in the classical background 
field. We will then use it in Eq. |3^ and the above equation to show that we do indeed 
recover the standard perturbative result for F 2 . For heavy quarks and/or moderate 
Q 2 , structure functions should be computed by inserting the fermion Green's func- 
tion in the definition for W^ v . We will see that in general the structure functions 
and the Fock space distribution functions are not the same. 

4 The Fermion Green's Function in the Classical Back- 
ground Field 

In this section we shall derive an expression for the fermion propagator in the clas- 
sical gluon background field described in section 2.2. The field strength carried 
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by these classical gluons is highly singular, being peaked about the source (corre- 
sponding to the parton current at x values larger than those in the field) localized 
at x~ = 0. Away from the source, the field strengths are zero and the gluon fields 
are pure gauges on both sides of x~ = 0. The fermion wavefunction is obtained by 
solving the Dirac equation in the background field on either side of the source and 
matching the solutions across the discontinuity at x~ = 0. Once the eigenfunctions 
are known, the fermion propagator can be constructed in the standard fashion. 
We begin this section with a discussion of the notation and conventions, proceed 
to write down the solution of the Dirac equation and finally, construct explicitly 
the fermion propagator in the classical gluon background field. This expression is 
formally exact and is valid to all orders in the source color charge density. 

4.1 The Dirac Equation in the Classical Background Field 

In order to compute the propagator for a spinor field in the fundamental represen- 
tation of the gauge group propagating in the background gauge field 

A% = 
A a _ = 

r-A t = 8(x-)K t (x t ), (37) 

where Kt(x t ),t = 1,2 is a two dimensional pure gauge, 

Kt(xt) = --V(x t )V t V\x t ), (38) 

we first need to solve the Dirac eigenvalue equation which can be written as 

[a t ■ (p t - gl t ) - V2p + a- - V2p~a + + /?m} ^ x {x) = Xtf; x {x) (39) 

for the spinor field ip and a corresponding equation for ip(x). The a's and (3 above 
are defined in appendix A and and p^ = —id^. For x~ < 0, the solution is trivial 

21 



and is just the free spinor plane wave solution. For x > the solution is less 



trivial and is given by the non-Abelian analogue of the 'Baltz' ansatz [39|. The 
full solution of the Dirac equation in the classical background field is 

= 0(-x-)e^ul x (q)+9(x-)^=J |^ dz t (V (x t )V\z t ))^ e^-^ x " 

{p 2 + M 2_ X) 



x e z H?t-p t ) exp 



2p- 



(«)■ 
(40) 



Above, the superscripts a, (5 denote the color index in the fundamental representa- 
tion and s the spinor index. The elementary spinors are normalized as u St \(q)u s ' t \i (q) 
2M5 ss i5\\t and summed over spins (uu)^ = (M — q 1 )^. 

The interested reader will notice that the above equation is not continuous 
across x~ = 0. This is because though ip-(x) is continuous across x~ = 0, V- is 
related to ^+ via the light cone constraint equation 



1 



V2q- 



a t -{-d t -gA t )+f3M 
% 



(41) 



which is discontinuous across x = by the same amount as in the previous equa- 
tion. 



4.2 Computation of the Fermion Propagator 

Having obtained the eigenfunctions for the Dirac equation in the classical back- 
ground field we are now in a position to compute the fermion propagator in the 
classical background field. This is given by the relation 

*<*■'> = /(f ^ + ^- te g*<*>ftM- <42) 

after identifying q + = (q 2 + M 2 — \)/2q~. It is straightforward to check from 
the above expression that (q 1 + M)S(x,y) = (2vr) 4 <5( 4 )(x - y). Substituting the 
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eigenfunctions from Eq. 40 in the above, we have 



S(x,y) = 8(-x-)e(-y-)So(x-y) + 9(x-)9(y-)(v(x t )So(x-y)VHy t )) 
A 1 ia-(x-v) f d2 Pt 



+ 



(2vr) 4 q 2 + M 2 — ie 



x e iPt-(x t -zt) el 



>q- 



(2n) 2 
(y{x t )V\z t \ 



d z z t d(x-)9(-y-) 



-{M 



+ 0(-x-)0(y-)e- ipt - {yt - Zt) 



exp 



. (pt + qt) 



2q- 



2q- 

■A- 



v(z t )vHyt 



x — (M-g/) 7 -(M-^-pi) 



where the free fermion Green's function is 

d A q iq . {x - v) (M-q 1 ) 



Sq(x - y) 



(43) 



(44) 



(2tt) 4 " q 2 + M 2 -ie' 

The translational symmetry of the Green's function in the x~ direction is of course 
broken by the presence of the source at x~ = 0. In the absence of the source of 
color charge, it may be confirmed that the free fermion propagator is recovered by 
putting V = I, where I denotes the unit matrix in the fundamental representation. 

The reader will note that the propagator between two points on the same side 
of the source, for either x - , y~ < or y~ > is the free propagator or a gauge 
transform of it. The only non-trivial contribution comes from the pieces connecting 
points on the opposite sides of the source. The 6{x~)0(—y~) piece can be written 
more simply as 

S(x,y) = -i J d 4 zV(x t )S (x,z) 7 -5(z-)S (z,y)V\z t ). (45) 

An analogous expression holds for the 9{—x~)9{y~) piece of the propagator. A 
similar simple expression for the scalar quark propagator |^(| was found recently 
by Hebecker and Weigert |^7| (see also the recent work of Balitskii |4^]). The only 
difference between the form of the above result and that for scalar quarks is the 7~ 
matrix present here due to the different spinor structure and a partial derivative 



d z + absent here due to the l/2g factor in Eq. 43 
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If we define 



G{x t ,x-) = 9(-x~) + 9(x~)V(x t ) (46) 

which is the gauge transformation matrix which transforms the gluon field at hand 
to a singular field which has only a plus component, A'^ = 5^ + a(xt), we then see 
that our propagator has the form 

S A (x,y) = G(x)So(x-y)tf(y)-iJd 4 zG(x)^e(x-)e(-y-)(V\z t )-l)- 

9{-x~)6(y-)(V(z t ) - 1) Jg%)S (s - z)-yS(z-)S (z - y) . (47) 

This very simple form of the propagator is useful in the manipulations below. 

In fact the current-current correlation function is explicitly gauge invariant. 
We may therefore use the singular gauge form of the propagator for computing the 
current-current correlation function 

St 9 (x,y) = Sv{x- y )-ijd 4 z\e(x-)6(-y-)(V\z t )-l)- 

0(-x-)9(y-)(V(zt) - 1) js (z " z)^ 5( Z -)S (z - y) (48) 

A diagrammatic representation of the form of the propagator above is shown in 
Fig. 1 In the expressions below for W^ v we will drop the superscript sing and 
simply use the singular gauge expression for the propagator. 

The Fock space distribution function however is gauge dependent. In computing 
it we must therefore either use the form of the propagator with the explicit gauge 
matrices above, or go back to our original form. In what follows, we shall use the 



original form of Eq. |43| for computing the Fock space distribution function. 

Our result for the fermion propagator in the classical background field was ob- 
tained for a ^-function source in the x~ direction. This assumption was motivated 
by the observation that small x modes with wavelengths greater than 1 / P + perceive 
a source which is a <5-function in x~ . The propagator above can also be derived for 
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Figure 1: Diagrammatic representation of the propagator in Eq. 48. 



the general case where the source has a dependence on x~ . The gauge transforms 
above are transformed from V(xt) — ► V(xt, x~), to path ordered exponentials, where 
V(xt, x~) is given by Eq. [2l|. Our result for the propagator is obtained as a smooth 
limit of Ax~ = l/xP + » x~(= 1/P + ). In other words, our form for the propaga- 
tor is the correct one as long as we interpret the ^-functions and 5-functions in x~ 
to be so only for distances of interest greater than 1/P + , the scale of the classical 
source. 



5 The Leading Twist Computation of 



Now that we have computed the fermion propagator in the classical background 
field, we are in a position to calculate the sea quark distributions in this background 
field, and in turn the structure function F2 and evolution equations using Eq. 
This calculation is accurate to lowest order in as but to all orders in asfJ,- Due to the 
singular nature of the propagators in the background field, the actual computation 
of the distribution function is a little subtle^and will be outlined below. 



25 



Before we go ahead to the computation, we will begin with a discussion of the 
averaging procedure over the labels of color charges at rapidities higher than those 
of interest. We will obtain a compact expression for it below. 



5.1 Color Averaging Over the Sources of Color Charge at Higher 
Rapidities 

Our expression for the propagator in the previous section makes no particular as- 
sumption about the color averaging over the color labels of the external sources 
corresponding to the valence quarks and/or gluons at higher rapidities than the 
rapidity of interest. The expression we quote before averaging is the quantity which 
will be useful in loop graph computations. To relate our expression to physical 
observables, as discussed in section 2 (see Eq. ||) we need to average over all the 
color labels of the external color charge density corresponding to the color charge 
density p a (xt,y) at rapidities greater than the rapidity of interest. 

Here we have smeared out the source in x~ and are no longer treating it as a 
delta function. This means that the sources acquire a rapidity dependence and that 
the weight for the Gaussian fluctuations over sources is replaced by 



Here the lower limit would be the rapidity of interest for evaluating the structure 
functions. 



labels corresponding to the partons at higher rapidities, we can employ the following 




(49) 



which leads us to define 




(50) 



If we average the Green's function in Eq. |43| over all possible values of the color 
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definitions for future reference. Denning 

^ < Tr (y(x t )V\y t )) > p = 7 (x t - y t ) , (51) 

we see that 

7(0) = 1, (52) 

which follows from the unitarity of the matrices V. Now defining the Fourier trans- 
form ^ 



7(pt) = J d z x t e-™* h(x t ) - 1] , (53) 

we have the sum rule 

The function "f{pt) will appear frequently in our future discussions and as we shall 
see, can be related to the gluon density at small x. 



For the particular case of a large nucleus [16] the averaging procedure has the 
form 



<0> p -- 



I»^(-f*/^)' (55) 



where fj? is the average color charge density per unit transverse area per unit ra- 
pidity. For an extensive discussion of the above averaging procedure we refer the 



reader to Refs. [16, 22, IS]. In appendix B we will explicitly derive an expression 
for j(pt) for Gaussian fluctuations. A Gaussian form for the averaging over color 
charges at higher rapidities is likely valid for very large nuclei or for realistic nuclei 
and hadrons for x « 1 but not at too small x's (or large enough parton densities) 
where non-linear corrections to the renormalization group equations are important. 



3 We define the Fourier transform in this way because it corresponds to only the connected pieces 
in the correlator. 
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5.2 Sea Quark Fock Space Distribution Function 



The relation Eq. [35] can now be combined with our expression in Eq. 43 to compute 
the sea quark distribution function. We shall use below the following identities: 



Tr[ 7 +(M-</)] = 4g H 



2q- 



: Tr [(M-c/-pl) 7 -(M-4) 7 + ] 



— (M 2 + q 2 +p t -q t ) 



(56) 



We then obtain the following expression for the distribution function 



dNf e 



2iN r 



d?k 



(fx dry 



drq e^i-^-^-v) 



(2vr) 3 J a J (2vr) 4 q 2 + M 2 — ie 

Uq+ [6{-x-)9{-y-) + 9{x-)9{y-) 1 {x t - y t )} 



x exp 

( Xpt + qt? -g 2 t _\ , , 

6XP V 2cp V p( Zt ~ Vt > 



9(x-)8(-y-)e ipt < Xt - Zt) 
j(x t - z t ) + e{-x~)6(+y~) e-^-te-*) 



(57) 



We will now sketch below the procedure used to simplify the above equation. 

a) First perform the integrals over the transverse co-ordinates. This introduces 
a common factor irR 2 and for the 9(—x~)9(—y~) term a factor 5^ (qt — kt) (a factor 



7 which was defined in Eq. 53 pops up in the other three pieces). 



b) Perform the integrals over x and y . For the 9(±x )9(±y ) pieces, we 
obtain the factors 



(58) 



q + — k + ± ie q + — k' + =F ' 

respectively. We have introduced above (to ensure smooth convergence) a slight 
difference in the momenta (k + and k' + respectively) multiplying x~ and y~ in the 



phases. In the final step we take the limit k 



k+ -> 0. 
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c) The simple contour integral over q + is done next. This introduces the factors 
0(±q~). 

d) The final step is to perform the (logarithmic) integral over q~ . The ultraviolet 
cutoff A — > oo cancels among the different terms and we obtain a finite result. 
Putting all the terms together and using the identity in Eq. [54] , we obtain the 
general result for the sea quark distribution function 

f dNf erm N c 1 r d 2 Pt 

ttR 2 dk+d 2 k t ~ 2^k+J (2^)2 l{Pt) 

' (k 2 + kf Pt + M 2 ) ( {k t + p t ) 2 + M 2 \ 
P 2 t+2k t .p t LOg { k 2 + M 2 ) ■ 

(59) 

In the region where the logarithm can be expanded, it can be checked analyti- 
cally that the argument of the above expression is positive definite. We have checked 
numerically that the argument remains positive definite in the entire (kt,pt) phase 
space (as it should be). 

In the next section we will study the above result in different limits and relate 
it to the well known evolution equations for sea quark distributions. 

5.3 Evolution Equations for Sea Quark Distributions at Small x 



In this section, we will show that for large q 2 , the Fock space seaquark distribu- 



tion we derived above in Eq. |59| gives us the Altarelli-Parisi evolution equation for 
seaquark evolution at small x. 

Towards that end, consider the Fock space distribution in the limit of large kt- 
In the integral in Eq. [59] we approximate kt » Pt- Then expanding the logarithm 



and defining 

(v.. _l ^2 _ j 

' -1 + K, (60) 



(k t + Pt ) 2 -k 2 
kf + m 2 
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we find that the terms in the square brackets [• • •] in Eq. ^ can be approximated 
by 

2 h-Pt , (h -Pt)K k 2 (k t -p t ) 



K K' 

2~T k 2 + M 2 ' k 2 + M 2 T k 2 + M 2 



+ 



(61) 



Now y(pt) has rotational symmetry in the transverse plane. This helps simplify our 
expression above since only even terms in kt • Pt survive. To leading order in p 2 jk 2 
then, Eq. |1 reduces to 



l dNf erm _ N c r d 2 Pt l 

irR 2 dk+d 2 k t 2W (2vr)2 7l;P ' j 2 



Pt 



k 2 t p 2 t 



ktP 2 t 



+ 



k 2 + M 2 3(k 2 + M 2 ) 2 {k 2 + M 2 ) 2 



(62) 



Since we are interested in the limit k 2 » M 2 , the above expression can be further 



simplified to read 



1 dNt< 



N r 2 1 



d 2 p t , . 



Pt 7(Pt) • 



(63) 
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To make contact with the evolution equations we will now obtain a relation 
between j(pt) and the gluon distribution function at small x. We begin with the 
relation we defined in the last section-Eq. p3|: 



lipt 



Then 



ph(pt 



d 2 x t e-^i^xt)-!}. 



d 2 x t ^d 2 xtl {x t ) 



(64) 



Recall that j(xt) = jj- < Tt(V(xt)V(0)^ > p . Expanding out the matrix V = 
1 + iA(xt) — A 2 (x t )/2 + •••", and doing the same for the correlator of gauge fields 
A 1 = ^-Vd l V^ , we obtain the relation 

in 1 



KJ 



PtliPt) 



2N r 



d 2 x t e ipt xt < Af(x t )A?(0) > p 



(65) 
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The correlator < AfAf > can be related to the gluon distribution function by 

the formula 

dN _ 2 | 1+ | 
~dfl ~ ~(2np 

Integrating both sides over l + , we obtain 



J d 3 xd 3 x'e-' a+x e +il+x " 9{x-)6{x'-)e ilt - Xt < A?(x t )A?(0) > p .(66) 



/ 



dx t e {A{x t )A{0)) - ^ R2 2 J ll+ld2l dl+ - 



Substituting the RHS of the above equation in Eq. |65| and the resulting expression 
for p 2 j(pt) into Eq. 63, we obtain 



dNf erm a s 2 1 f 1 dy f k t , , dN9 lue . . 

x 1.11.2 = T7-oT2 / — / dPtX-r-rr, (68) 



dxdk 2 2ir 3 k 2 J x y Jo dydp 2 
where we defined x = k + /P + and y = l + /P + . 

Recall that the structure function to leading twist and lowest order in as is 

dN ferm 

/o dxdKi 



F 2 {x,Q 2 )= dk 2 x = x[q{x,Q 2 )+q(x,Q 2 )} = 2xq{x,Q 2 ), (69) 

Jo dxdkf 

and similarly, 



xG(x,Q 2 )=[ x—j, (70) 



lo dxdkf 

where xq(x,Q 2 ) and xG(x,Q 2 ) are the quark and gluon momentum distributions 

respectively. At small x, the sea is symmetric and we take q(x,Q 2 ) = q(x,Q 2 ). In 

the limit of light quark masses Q 2 >> M 2 we find then that 

d(xg(x,Q 2 )) _as2 f 1 dy q2 , m s 

dlog(Q2) ~ An3 X J x tfVGW'Q)' ^ 

Now at small x, yG(y, Q 2 ) is slowly varying. One can for instance parametrize it by 

a power law 6v» s > where C is some constant. The scale of variation of the structure 

function then corresponds to higher orders in as- We can therefore take yG(y, Q 2 ) 

out of the integrand. (The same is true for any other slowly varying function |36|.) 

We then get finally for the seaquark evolution equation at small x the result 

d{xq(x,Q 2 )) a s 2 2 

W^3 lGM)+ ° (as) ' (72) 
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Thus to lowest order in as, the seaquark evolution at small x is local and simply 
proportional to the gluon density at that x. 

Now consider the Altarelli-Parisi evolution equation 



d£ _ a s (Q 2 



dQ 2 2vr 
Above the operation ® denotes 

A® B = 



[V®P qq + G® 2fP qG ] 



1 dy . , , „ , x , 
-A(y)B(- 

y y 



(73) 



and £ = ^2f(q + q), where / is the number of flavors. Also, P qq and P q c are the well 
known Altarelli-Parisi splitting functions. Since at small x the quark distribution is 
as suppressed relative to the glue, taking q = q and / = 1 the leading contribution 
to the seaquark evolution is 



d(xq(x,Q 2 )) 1 as f 1 dy 



(74) 



dlog(Q 2 ) 

Above we have made use of the relation 2fP q o(z) = f[z 2 + (1 — z) 2 ]. Let z = y/x. 



Then the above relation can be re-written as 



d(xq(x,Q 2 )) a s 



dz 



. „ zG(zx, Q ) 

47T Jl Z 2 



+ 1 



(75) 



d\og(Q 2 ) 

Again, as previously, we can argue that since at small x zG(zx, Q 2 ) is slowly varying, 
we can take it out of the integral. Doing that and performing the integral, we obtain 
finally 

d(xq(x,Q 2 )) a s 2 



-xG(x,Q 2 ) + 0(a s ). 



(76) 



dlog(Q 2 ) 4vr 3~ 

which is the same as Eq. fl%. This form of the sea quark evolution equation at small 



x was first obtained by Ellis, Kunzt and Levin [35|. 
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6 Computation of Current- Current Correlator to All 
Twists in the Classical Background Field at small x. 



In the previous section we used the fermion Green's function derived in section 
4.2 to compute the lightcone seaquark distributions at small x and subsequently 
the leading twist expression for the structure functions. It was shown that these 
structure functions obeyed evolution equations which were precisely the small x 
Altarelli-Parisi evolution equations. 

In this section we will again use the fermion Green's function in Eq. ^ to 
derive an explicit result for the hadronic tensor W^ u . This result will be valid 
to all twists at small x and for arbitrary quark masses. For light quarks, we will 
compute the structure functions F\ and F2 and show that the leading twist result 
in the previous section is recovered as a limit of our general result. We will also 
use our general result to obtain expressions for heavy quark structure functions at 
small x. We should note here that Levin and collaborators have studied screening 
corrections to the structure functions for light and heavy quarks in the Glauber- 
Gribov framework 



6.1 Analytic Result for at Small x 

As in the previous sections, we define 

W^ v (q,P,X-) = lm [ d 4 ze iq - z <T(J»(X- + -)J U (X- --)>, (77) 
J 2 2 

where "Imaginary" stands for the discontinuity in q~ . Then 

W>* u (q,P) = —a— [ dX-W^(q,P,X-) = —aP + lm [ dX~ [ d^ze 1 ^ 
2i ra J 2n J J 

' ^{S x A,X- + ^X-- Z -]ySaJx-- Z -,X- + Z -) ,") . (78) 
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The only terms in the propagator that contribute to the above are the 0(±x )0{ z py 



pieces. Using our representation for the propagator in Eq. 45 



W^(q, P, X) = lmN c J d*z i** J d\ t d^ t (2n)5(p- - l~) 

x (2ir)5(p'~ - i'-yv\x^-u t ) e _«.(x-f-ut) e i'-(x-i-u') t e - ip '.(x+i~u> t ) 

Tr f jM-fyrVj - f)r(M - l'/) 7 -(M-p'/)Y \ 

{{p 2 + M 2 - ie) {I 2 + M 2 — ie) (I' 2 + M 2 — ie) (p' 2 + M 2 — ie) J 

x (o{x- + z -)e( z --x-) + e(x-- z -)e(-x-- z ^ 1 {u t - u ' t ). (79) 

We have used above the condition ^{ut-,u' t ) = 7(u£ — Ut) = j(ut — u' t ), since the 
correlation functions are translationally, rotationally and parity invariant in the 
transverse plane. We now perform the integral over X~ and use the identity 

J dx- e n-^»-'+^)x- (e ( x- + |>»c§ -x-) + Hx- - f )«(- x- - 1) 

= <^ {P + - ,+ -V + n 2sin - i+ + 1 ' + ~ *' + >) ■ <80) 

We then obtain 

W^ q ,P) = ^-ImJd^*e(z )J—— F ^8(p -I )S(p 
x (2^) V 2 )(Z t - w + pj - /0 ^ { l l^HzA±A) e i-Mn+h- P ' t -i' t )/2 e i Z +(p-- P '-) 



x Tr 



(M - j/)7"(M - r)7 M (M - l'/h~(M - pfYf 
(p 2 + M 2 — ie) (I 2 + M 2 — is) {V 2 + M 2 — ie) (p l2 + M 2 — ie) 

1 2sm( Z -(p + -1+ + -p' + )) . 



(p+ - /+ -p /+ +/ /+ ) \ 2 
The subsequent procedure of solving the integrals is as follows: 

a) Perform first the integral over li. Then perform the integral over zt and z + . 
Defining k t = (p t - l t + p' t - l')/2, this sets l' t = p t - h - q t , h = Pt ~ h and 
qt = Pt — Pt- Also, we get q~ = p~ — p' . 

b) Next perform the integral over z~ . This gives us 

' dz-e(z-) e^'~ - e -^'d + -i' + )) 
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l' + -l+-p' + +p+ 

{q+ -p+ +p ,+ + ie){q+ -1+ + V + + ie) V ; 
The numerator above cancels the term l/(p + — l + — p' + + l' + ) in W^ v '. 



(82) 



c) Lastly, we do the integrals over p' + and l' + . This sets p' + = p + — q + and 
l' + = l + — q + . Then we can define in , p' = p — q and I' = I — q with I = p — k 
and k~ = 0. 

After these considerations, we can write W^ u as 

tlrliw „. aP + N CT f d 4 p d 2 k t dk + , . 

W^(q,P) = -Im / — Vj- — ~ — -7W 

v ; 2vrm ./ (2vr) 4 (2?r) 2 (2vr) ,v ; 



x Tr 



(M - p 1 )^- {M-IYf (M - I'/)*/- (M - p'/)Y 
(p 2 + M 2 — fej (Z 2 + M 2 — ie) (Z /2 + M 2 - ie) (p /2 + M 2 — ie 



(83) 



where I', p' and Z are defined as in step 'c' above. Correspondingly, we can write 
W^ v as the imaginary part of the diagram shown in Fig. 2. 

For the case of deep inelastic scattering, q 2 > (see footnote 1), and we can 
cut the above diagram only in the two ways shown in Fig. 3 (the diagram where 
both insertions from the external field are on the same side of the cut is forbidden 
by the kinematics). 

Also interestingly, the contribution to W^ v can be represented solely by the 
diagram in Fig. 4 and not, as is usually the case, from the sum of this diagram and 
the standard box diagram. This is because in our representation of the propagator 
multiple insertions from the external field on a quark line can be summarized into 



a single insertion. See for instance Eq. 48 which makes this point clear 



Applying the Landau-Cutkosky rule, and making the shift p — > p + k, Eq. 
can be written as 

W ' ' 27rm J (2tt) 4 (2tt) 2 (2tt) 71 * J 
x Tr { (M — p 1 — y)7 _ (M - rf)-f{M + 4)-y~(M - jf - kf + j)j u } 
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/ 

/k 




Figure 2: Polarization tensor with arbitrary number of insertions from the classical 
background field. The wavy lines are photon lines, the solid circle denotes the 
fermion look and the dashed lines are the insertions from the background field (see 
Fig. 1). The imaginary part of this diagram gives W^ v '. 



0( p + + k +)0(q + - p + )(2ir) 2 5({p + kf + M 2 ) «5((p - qf + M 2 )) 



p 2 + M 2 



(p+k-q) 2 + M 2 

1 1 



+ e(p + )6(q + - p + )(2ir) 2 5(p 2 + M 2 ) 5((p + k - q) 2 + M 2 )) 

(84) 



(p+k) 2 + M 2 {p-q) 2 + M 2 

With an appropriate change of variables, the second term is the same as the first 
except that now fi <-> v. We then get 

aP+N c f d 4 p d 2 k t dk+ 



W^(q,P) = 



2nm J (2vr) 4 (2vr) 2 (2tt) 



^{k t )M^e{p + + k + )6(-p + ) 



x (27r) 2 5((p + k) 2 + M 2 )5({p-q) 2 + M 2 )) 



p 2 + M 2 (p + k- q) 2 + M 2 ' 



(85) 
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/ k 



p-q 




\k 



Figure 3: Cut diagrams corresponding to the imaginary part of W^ v . 
where above the trace is represented by |^| 

M» v = Trj (M - p 1 - V)7" (M - jf)^{M - + (J)j- (M - j/ - k/ + 4)^ u + y. <-» 1/ 



Performing the (5-function integrations above, our expression for VF' t!/ can be sim- 
plified to 



(86) 



-a 



P+N r 1 



167T771 ((/ 



~V2 



(2vr)2 (2tt)2 7(fc * j 



2 9 - 



dp+ M^ v 



I(h,p t ,q,p + ), (87) 



with the definitions M 2 _ q 



(pt ~ Qt) 2 + M 2 and M* +h _ q = (fit + h~ Qtf + M 2 , 



4 Kinematic note: the observant reader will notice we have put q + — here. Since we are 
working in the infinite momentum frame, the hadron has only one large momentum component, 
P + . The rest are put to zero. For the photon, we choose a left moving frame such that q° = \q z \ 
and q + = 0. Then, q 2 = q 2 > 0, P ■ q = -P + q~ and x B j = -q 2 /(2P ■ q) = q 2 /(2P + q~). Since 
in the infinite momentum frame < xsj < 1, this gives q~ > 0. We are at liberty to choose the 
above frame since the hadron tensor is clearly Lorentz invariant and hence can be expressed purely 
as a function of q 2 and P ■ q. The explicit presence of P + in Eq. ^ may give the reader cause 
for concern. It arises from a relativistic normalization of the vacuum states. One can show that, 
despite appearances, Eq. [79 is Lorentz invariant. Of course our later results will confirm this fact. 
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Figure 4: In the singular gauge representation for the propagator (see Eq. 48 and 
Fig. 1), multiple, higher twist contributions from the classical gluon background 
field to the current-current correlator (imaginary part of LHS) is equivalent to the 
imaginary part of RHS. 



and 



1 



P 2 1~ P (29-)A^ +fc _, 



Eq. S7 is our general result for the hadronic tensor. We shall now study the different 
components of the above tensor and extract from these, different limits of interest 
to us. 



6.2 Structure Functions at Small x 



The hadronic tensor W^ u can be decomposed in terms of the structure functions 
Fi and F 2 as p 
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where is the four-momentum of the hadron or nucleus and P 2 = m 2 ~ 
(<< q 2 ). In the infinite momentum frame, we have P + — > oo and P~ , Pt ~ 0. 
The above equation can be inverted to obtain expressions for F\ and i*2 in terms of 
components of W^ v . Since in our kinematics g + = (see footnote 3 for a kinematic 
note) 

/ „2 \ 



with 



+ 17^5 ^"". ( 9 °) 



It is useful to verify explicitly that our expression for W^ v derived in an external 



field can be written in the form of Eqn. 89. Recall that W^ v can be written in 
Lorentz covariant form by using the vector = <5^ + . Using n-j = —7" in Eqn. 
we see that W^ v is a Lorentz covariant function of the only vectors in the problem- 



and n^. Identifying = P^/P + in Eqn. 89, we see that these forms are in 
complete agreement. We also see that all factors of m disappear from F\ and F2 



by the explict forms of Eqns. p9| and S3. Henceforth we will take m = 1 since it 



disappears from the quantites of interest, and was in fact only introduced due to 
historical normalization conventions. 

We also see that the structure functions can only be functions of q 2 and n ■ q 
by Lorentz invariance. We can therefore take q + = for the purpose of computing 
F x and F 2 . 

To compute iy ++ and W , we need to know the the traces M ++ and M , 



respectively in Eq. |87[ We can compute them explicitly and the results can be 
represented compactly as 

-3-M++ = \ (M 2 M 2 +k _ q + M 2 +k M 2 _ q - q 2 k 2 ) , (92) 
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and 



M— =32(p-) 2 (p- -q~) 2 . (93) 

From the relations above of F\ and F2 to and W , we obtain from Eq. ^ 

the following general results for the structure functions for arbitrary values of Q 2 , 
M 2 and the intrinsic scale /U, [] 

F 2 = ^k 1*****,^ Mil log ( M U M ^ \ m 

and 

Above, we defined the function J 7 in terms of a and /? (in turn functions of pt, kt, qt 
and M) which are defined as 

The general expression for is 



1 3 1 1 ^ d j_ i 2 

3 ~ 2 ( a + ^ + W + a~S + /? 



(.97) 



Eqn. ^ and Eqn. |95] are the central results of this work H. Since they are the most 
general possible expressions for arbitrary values of Q 2 , M 2 , and // 2 , it is inevitable 
that they look complicated. We shall show in the following sub-section that they 
simplify considerably in the high q 2 limit. 

5 which is implicitly contained in the function j(kt) in Eq. 

6 Our expression for F2 in Eq. can be further simplified-a result which will be discussed in 
a forthcoming paper E3]. 
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6.3 Structure functions in the limit q 2 — > oo. 



We shall now obtain the leading twist limits of Eq. 95 and Eq. 94. In particular 
we will show that our structure function for M 2 — ► and q 2 — ► oo is identical to 
the structure function obtained by integrating the light cone Fock distribution in 



Eq. 59. That this should be the case is a well known property of the leading twist 



structure functions [25, 26]. Further, we will recover the Callan-Gross result 
F\ = F2/2X in this limit. 



Consider first our general formula for F2 in Eq. 94. The trace simplifies consid- 
erably when we put M = 0. For qt 3> Pt,kt, we obtain 



In the logarithm, the ratio M 2 _ q /M 2 +k _ q 
integrals, 



16q 2 (p 2 + P fk t ) . 

1. Finally, in the denominator of the 



{mI_„mI 



Mi 



Mr 



M 2 (Mj 



2 

p+k 



AC 



*p-q ±yl p+k lvI p+k-q^"p 

Then putting these back into our general expression [] we obtain 

oN r r , 9 r d 2 k, 



2vr 4 



(2vr) ; 



1 _ (Pt+h-pt) lQg / (fct + p 4 ) 2 



(A; t 2 + 2fe t • ft ) 



(98) 



A comparison with Eq. [5^ immediately reveals that, setting kt Pu an d integrating 

the latter over pt upto q 2 gives an identical result to the one above. Thus we have 

recovered a well known, non-trivial leading twist result as the limit of our general 

expression for q 2 — > 00 and M — > 0. 

Even though our general expression for the longitudinal structure function F\ 

looked terribly complicated, in the limit considered here it is remarkably simple. 

In this limit a, (3 — > 1 and hence Eq. ^ for T — * 3— a constant! The product in 
7 The contribution in the pt integral in the region p t ~ qt is identical to that from pt << q t . 
This provides a factor of 2 that must be taken into account. 
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the denominator of Eq. 95, M 2 _ q M 2 +k _ q — ► q 4 and cancels the q 4 factor outside 
the integral. From the sum rule Eq. 54, we find remarkably that the complicated 



integral vanishes and Eq. 95 reduces to 

fl = g. (99) 

The above is the well known Callan-Gross relation. 

We should clarify the result obtained above to avoid confusion. The reader may 
note above that the deviation from the Callan-Gross relation vanishes as a power law 
as q 2 — ► oo. On the other hand, it is well known in QCD |28|, 29 1 that the violations of 
the Callan-Gross relation only disappear logarithmically as q 2 — > oo. The apparent 
contradiction is resolved by one realizing that the logarithmic violations at large q 2 
in QCD come from diagrams where the sea quark emits a gluon (thereby violating 
Feynman's parton model helicity argument). These diagrams are of higher order in 
our picture and are therefore not included. In fact, the deviations from the Callan- 
Gross relation of the sort discussed above (at small x) should die off faster than 
logarithmically at very large q 2 because for sufficiently large q 2 , the violations of 
the Callan-Gross relation should come from precisely the diagrams not included 
here. At moderate q 2 however, the contributions we have discussed above should 
be important. 



7 Summary 



In this paper we have used a classical theory of the gluon field to derive expressions 
for Fock space distribution functions of quarks and structure functions for deep 
inelastic scattering. This theory is valid at small x when the gluon density is large. 
In this region, the coupling constant evaluated at this density scale is small. With 
this density scale denoted by fi 2 , we have seen that when q 2 » a 2 fx 2 and q 2 » M 2 , 
where M is the mass of the heavy quark being probed, all leading twist results are 
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reproduced. For q 2 < a 2 fi 2 or q 2 < M 2 , we derived an expression valid to all orders 
in twist but only to leading order in as- In this kinematic region, the Callan- 
Gross relation is not valid, and there is no simple relation between the Fock space 
distribution function and the structure functions for quarks. 

The structure function of heavy quarks deserves more study. If M 2 >> a 2 fi 2 , 
then for q 2 » a 2 fjL 2 , the heavy quark distribution is a linear function of the gluon 
distribution function. The gluon distribution function may nevertheless be com- 
puted to all orders in twist. The leading twist relation between the structure func- 
tion and the quark distribution function is however not valid for large M 2 . As we 
go to smaller values of x corresponding to larger values of /i 2 , the non-linearities 
corresponding to higher powers of the gluon density turn on and can be studied 
systematically in our weak coupling formalism. 

The situation for light quarks is also amusing and needs more study as well. 
If q 2 < a 2 [i 2 , then the non-linearities in the gluon distribution function become 
important. In this kinematic region, we expect saturation of the gluon distribution 
function, and our function "y(kt) ~ up to logarithms of kt- If this is the 

case, a look at the definition of F2 and F± shows that these distributions should 
be dimensionally of order q 2 . In this region a precise analytic estimate is difficult 



since in the integral representations for F\ and F2 (Eqs. |95| and 94 respectively), 
the there is no hierarchy of momentum scales. All momenta are of order q, and the 
integrand does not simplify much. Nevertheless, we see that the saturation of the 
gluon density is sufficient to imply saturation of the quark density. 

Both the study of the heavy quark and light quark distributions merit more 
theoretical and phenomenological work within the framework described in this pa- 
per. 
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Appendix A: Notation and Conventions 



We start by defining our convention and notations. Our metric is the +2 metric 



(—,+,+,+). The gamma matrices in space-time co-ordinates are denoted 



by carets. In the chiral representation, 



7 



/ 



V 



o / 

/ 



; f 



/ 



V 



a 1 
-a* 



7 



/ 
—I 



and {7^,7 l/ } = —2g^ u . Above, a l ,i = 1,2,3 are the usual 2x2 Pauli matrices and 
/ is the 2x2 identity matrix. In light cone co-ordinates, 7 =t = (7 db 7 3 )/\/2 and 
{7^,7^} = —2g^ u , where g ++ = g = 0, = g h = —1 and gt x t 2 = 1 where 
tx,t2 = 1,2 here stand for the two transverse co-ordinates. Note for instance that 
in this convention A + = —A~ and At = +A t . Also, q 2 = —2q~q + + q 2 hence a 
"space-like" q 2 implying large space-like components would correspond to q 2 > 0. 
We now define the projection operators 



a 



V2 



(100) 



which project out the two component spinors Some relevant 
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properties of the projection operators are 



(o^f = a ± ; o^cF = ; a + + a~ = 1 ; (a ± ) 



(101) 



It follows from the above that if) + + ?/>_ = ^. In the following, we will also use the 
familiar Dirac conventions (3 = 7 and a± = 7°7^. 

The two component spinor is the dynamical spinor in the light cone QCD Hamil- 
tonian Pq CD and it is defined in terms of creation and annihilation operators as 

d 3 k 



fc+>0 2V4( 27r )S 



[e ik ' x w(s)b s (k) +e- ik - x w(-s)dl(k)] . (102) 



Above b s (k) is a quark destruction operator and destroys a quark with momentum 
k while d\(k) is an anti-quark creation operator and creates an anti-quark with 
momentum k. Also above the unit spinors w(s) are defined as 



«4> 



/ \ 
1 


v y 



/o\ 



1 



(103) 



The creation and annihilation operators obey the equal light cone time (x + ) com- 
mutation relations 

{b s (k, x + ), b\, (k>, x+)} = {d s (k, x+), 4 x+)} = (2tt) 3 <5( 3 ) (fc - fc')<W . (104) 

The above definitions ensure that the light cone QCD Hamiltonian can be defined as 
Pq CD = Pq + Vqc d , where the non-interacting piece of the Hamiltonian is defined 

as 

tr'k ,-, (fc| + M 2 ^ 



P ~ J ( t)3 E 



2fe+ 



(&^M*0+4(*K(*0) 



(105) 



The definition of quark distribution functions is further discussed in the text of 
section 3.2 
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The dynamical components of the gauge fields Af(x) with i = 1, 2 in light cone 
gauge A + = are defined as 

A?(x) = / £ ^ e tt -oS(fc) + e- ife -af (fc) , (106) 

where the A's here correspond to the two independent polarizations and f (a") 
creates (destroys) a gluon with momentum k. They obey the commutation relations 

[a%(k), a b y \k')] = (27rf5^(k - k')5 ab 5 xy . (107) 



The gluon distribution function is then defined as 

dN _ af a? 
d?k ~ (2tt) 3 ' 



(108) 



Performing the Fourier transform of the gauge field above, we obtain Eq. 6(: in 
section 5.3 For a more extensive discussion of the above formalism see Ref. f42| . 

We should mention here that there are several conventions in use. For a discus- 
sion of some of these, see the review article by Brodsky, Pauli and Pinsky [43|. Our 



convention is alike that of Kogut and Soper |40j but differs from theirs for quark 



spinor and gauge field normalizations by a factor vT^+T/ \J (2tt) 3 . 

Appendix B: Derivation of the Function ^(pt) for Gaus- 
sian Fluctuations 

Since j(pt) is defined as the Fourier transform of ^ < Tr(U(xt)U^ (yt)) > p , we 
need to compute this correlator in co-ordinate space first. Note that the symbol 
<•■■>/> denotes the averaging over with a Gaussian weight. Now, as discussed by 
Jalilian-Marian et al., if 



U(y,x t ) = U ooy {x t ) = Pexp 



i / dy'A(y',x t ) 



ii 



(109) 
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where U OC: y(x t ) is the path ordered exponential (in rapidity) which corresponds 
to the pure gauge potential A 1 = —U 00 ^{xt)V L Ul Q ^ y {xt)/ig, then A % satisfies the 
Yang-Mills equation 

j Ai,a 

D l -—=gp a (x t ,y), (110) 
dy 

if A, the argument of the path ordered exponential, satisfies the Laplace equation 

V 2 A a (x t ,y)=p a (x t ,y). (Ill) 
The measure for the functional integral is then 

m exP (- jf" * / - fw exp (- f * / **^^) m 

As argued by Jalilian-Marian et al., we can write 

EWxt) =: U^xt) : exp ^ 4 N c r(0) ; (n3) 

where : • • • : denotes normal ordering and 
1 x 2 

r(x t ) = ^4 = T(0) + ^ log(x 2 L 2 ) + finite pieces... . (114) 

Above, r(0) cx 1/L 2 where L corresponds to an infrared cut-off. Writing Uoo, y {xt) in 
the above normal ordered form enables us to isolate and exponentiate the infrared 
singular terms coming from disconnected graphs. 
Our correlator has then the form 

j{x t ,x t ;y,y) = N c j'(x t , x t ; y, y) exp (^- 9 Nc ^^ + j^j dy'u 2 ^')^ ,(115) 

where 

r2 Tr(AV 4 A)' 



7 



// f°° f TriAV A) \ 
[dA]exp {- J Q d V J d 2 x t g4 ^ 2 ^ J ( : U °°,y( x t) ( : ^oo,y(zt))' • 

(116) 
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Expanding out first few terms in the path ordered exponentials above, we have 



/(i) 



1, 

N 2 — 1 f°° 

-T(x t - x t )g A 9{y - y) dy'^ 2 {y'). 



2N r 



1 



/(2) 



1 

2! 



/V 2 — 1 f°° 



1 2 



(117) 



In the expression for j'^ above only one of the two possible terms survive on 
account of the path ordering. From similar considerations it can be argued that in 
general 



2N, 



T(x t - x t )g 4 £(y) 



(118) 



where £(y) = fjj° dy'[i 2 (y'). Resumming the terms above and including the discon- 
nected pieces, we have 

7 (x f - x t ; y, y) = exp ^ (T(x t - x t ) - T(O))^ . (119) 

The above expression is the complete non-perturbative result for ~f{x t — xt;y,y')- 
Note that trivially 7(0; y, y) = 1 as we would expect from the definition of 7. 
Taking the Fourier transform of 7, 



lipt 



d 2 x t e ipt - Xt j{x t ) 

d 2 x t e ipt - xt exp [K(T(x t ) - T(0))] 



(120) 



where k = g {N* — l)£(y)/2iV c , and expanding out j(xt), we obtain 



5 {2) (pt) + kJ d 2 x t e ipt - x *T(x t ) - KT(0)6^( Pt ) + ■■■ 



(121) 



If we now recall the definition of T(xt) from Eq. 114 

d 2 k t e ikvxt 



r(x t ) 



(2vr) 2 kf ' 



(122) 
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and substitute for T(xt) in the above, we find 

■m = (i - *w**<*d + ^ ( -^ + - . (123) 

The first term in j(pt) is not relevant for our computation of the sea quark dis- 
tributions since the relevant momenta are pt » Aqcd- The second term is the 
perturbative expression computed by us previously ^(|-upto a factor of two which 
was missing in that paper. Note also that the perturbative ^(pt) above explicitly 



satisfies the sum rule condition of Eq. 54. 
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